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Abstract. 

We present the results of a search for a high-energy axion emission signal from ^Li 
(0.478 MeV) and D(p,7)'^IIe (5.5 MeV) nuclear transitions using a low-background 
7-ray calorimeter during Phase I of the CAST experiment. These so-called "hadronic 
axions" could provide a solution to the long-standing strong- CP problem and can be 
emitted from the solar core from nuclear Ml transitions. This is the first such search 
for high-energy pseudoscalar bosons with couplings to nucleons conducted using a 
helioscope approach. No excess signal above background was found. 

Keywords: axions, axion-photon coupling, axion-nucleon coupling, hadronic axions 
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1. Introduction 

The observed CP invariance in the strong interactions is not a priori expected, as 't 
Hooft pointed out [1], and has been named the strong-CP problem. Nonperturbative 
effects in the theory give rise to a CP violating "6'" term which appears in the QCD 
lagrangian as 

Ce = e^C^^G'-". (1) 
ovr 

Here, G^j, is the gluon field strength and G^'^ = |e^i,po-G''°" its dual. The apparent CP 
invariance of QCD derives from the fact that 6 is measured to be vanishingly small via 
the neutron electric dipole moment, for which the current upper limit is \dn\ < 6.3x 10~^^ 
e cm [2]. This limit on dn implies an upper limit 6 < 10~^° [3]. 

In 1977, Peccei and Quinn proposed a physical origin for ^ = by introducing a 
global U{1) chiral symmetry [4], often referred to as f/(l)pQ. The parameter 6 thus 
becomes a dynamical variable that is forced to zero when the potential is minimized. 
Weinberg and Wilczek showed that such a solution implies the existence of a new 
particle, the axion, and that such a particle can have couplings to quarks, nucleons, 
leptons and photons [5l [6] . 

The axion was first thought to have couplings on the order of the weak scale [5] and 
a mass of ~200 keV. Experimental evidences against this coupling strength and mass 
range, most notably through limits on the magnetic moment of the muon, kaon decay 
and quarkonium studies, prompted the idea of "invisible" axions. There are two classes 
of invisible axion models: KSVZ (Kim, Shifman, Vainshtein, and Zakharov) [71 E] and 
DFSZ (Dine, Fischler, Srednicki, and Zhitnitskii) [9l HO]. In the former axion models, 
commonly referred to as hadronic axion models, couplings to leptons are strongly 
suppressed. Since couplings to nucleons and photons remain, detection of the hadronic 
axion is still possible. Other models have also been proposed with suppressed axion- 
photon coupling [11], but we focus here on axion models which include couplings to 
both photons and nucleons. 

Most of the experiments searching for axions or similar pseudoscalar bosons [12] 
have been relying on the coupling to two photons i.e. Primakoff effect [I3]. The axion- 
photon coupling is given by the effective Lagrangian 

C,^ = -^g.^F^'''F^,a = g,^E-Ba, (2) 

where a is the axion field, F^'^ the electromagnetic field strength tensor, and F^^ its 
dual, E the electric and B the magnetic field of the coupling photons. The effective 
axion-photon coupling constant g^-y is given by 



a 



2vr/a 



E 2(4: + z + w) 



'--1.95 ±0.08 , (3) 



_N 3(l + ^ + w; 

where /a is the Peccei-Quinn symmetry breaking scale, z = m^/m^ = 0.56 and 
w = m^/nis = 0.028 are the quark mass ratios. Here E and are the model dependent 
coefficients of the electromagnetic and color anomaly of the axial current associated with 
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the ?7(l)pQ symmetry, respectively. Frequently cited axion models use E/N = [TJ |B] 
or E/N = 8/3 [9l |T0] but, in general, E/N can take different values depending on the 
specific model details. 

In the Primakoff process an axion couples to a virtual photon in an electromagnetic 
field and converts to a real photon, or vice- versa. The conversion probability in vacuum, 
in a magnetic field of length L and strength B, depends on the axion-photon coupling 
constant g^^r^ and the momentum transfer between the axion and the photon q as given 



In addition, this probability depends implicitly on the axion mass and the axion 
energy E'a through the relation q = m'^/2Ea. For gL < 1 the conversion is coherent and 
equation (j4]) reduces to {gs^^BL/2Y . 

Axion-photon mixing permits a variety of production mechanisms and detection 
techniques, many of which were first pointed out by Sikivie in 1983 [15]. Magnetically 
induced vacuum birefringence [TB], stellar and terrestrial magnetic fields, pulsar 
magnetic fields, and resonant cavities [15] all provide methods for the production or 
detection of axions. Experimental and astrophysical limits on the axion are typically 
stated in terms of the photon coupling, g^^, versus mass, ma, and recent experimental 
and cosmological limits are shown in figure [H 

Due to the axion coupling to nucleons, there are additional components of solar 
axions emitted in nuclear de-excitations and reactions. The energy of these mono- 
energetic axions corresponds to the energy of the particular process. In this paper 
we present the results of our search for mono-energetic solar axions which may be 
emitted from ^Li* de-excitation and D(p, 7)^He reaction by using the CERN Axion 
Solar Telescope (CAST) setup. To detect photons coming from conversion of these 
axions in the CAST magnet, we used high-energy photon calorimeter that was mounted 
on one end of the magnet during the CAST phase I. 

In section |2] nuclear axion couplings are discussed and the signal expected from 
nuclear axion emission in the Sun from the first excited level of ^Li (0.478 MeV) and 
reaction D(p, 7)^He (5.5 MeV) is described. In sections [3] and H] the data selection, 
systematics and analysis of the data are then presented in detail. 

2. Nuclear axion emission in hadronic axion models 

2.1. Axion-nucleon coupling 

Coupling to nucleons occurs through the spin operator a [26] and because axions 
carry spin-parity = 0", 1"^, 2~, ... nuclear deexcitation via axion emission occurs 
predominantly via Ml magnetic nuclear transitions. Several channels exist for solar 
axion emission via these transitions [271 EHl EHl [301 ED 132] • Here we focused our attention 
on the thermonuclear fusion reaction p + d— j-'^He + a (5.5 MeV) and associated reaction 



in [H] 




(4) 
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Figure 1. Exclusion plots in the axion-photon coupling versus the axion mass plane. 
The limit achieved on PrimakofF axions from the X-ray detectors by the previous 
Phase I of the CAST experiment (updated limits consistent with expectations are 
reported in [18]) is compared with other constraints (Lazarus et al. [19], SOL AX |20j . 
COSME [21], DAMA [22], Tokyo helioscope experiment SUMICO [23] and HB 
stars [24]). The vertical red line (HDM) is the hot dark matter limit for hadronic 
axions ma < 1-05 eV [25] inferred from observations of the cosmological large-scale 
structure. The yellow band represents typical theoretical models with \E/N — 1.95| in 
the range 0.07-7 while the green solid line corresponds to the case when E/N = is 
assumed, as in the KSVZ model [7l[8]. 



chain ^He(a, 7)^Be(e~, i/c)^Li*— j-^Li-l-a (0.478 MeV) as a source of solar mono-energetic 
axions. 

The branching ratio for axion emission via Ml transitions is directly calculated 
as [26] 

ra_ 1 (k.Vf 1 \ 

Here, kg, and are the axion and photon momenta, respectively, and will both be 
approximately equal to the decay channel energy. In addition, a = 1/137, 6 is the 
E2/M1 mixing ratio, /iq ~ 0.88 and fis ~ 4.71 are the isoscalar and isovector nuclear 
magnetic moments, respectively, in nuclear magnetons. 

The remaining terms r] and /3 are nuclear structure dependent parameters which 
depend directly on the initial and final state nuclear wavefunctions. It is via these terms 



go(3 + 93 
(/io - l/2)/3 + /is - r/ 



(5) 
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that the decay- channel specific axion-physics is expressed, and are given by 



V 



(J/IIEf=iiWr3WI|J.) 



(3 



|J.) 



where Jj and J/ are the angular momenta of initial and final state respectively, 

is the nucleon orbital angular momentum operator, a{i) represents the spin operator, 

while T3{i) is the isospin operator. 

In equation (jS]) we have written the axion- nucleon coupling as {gofS + gs), where go 
is the isoscalar coupling and g^ the isovector coupling. In the hadronic axion models 
these are written as [H] 

l + z-2w' 



90 



93= - 



mN 1 

mN 1 
X2 



2S 



(D + F) 



(3F - D) 
1 - 



1 + z + w 



1 + z + w 



(6) 



where S = 0.4 characterizes the flavour singlet coupling, F = 0.460 and D = 0.806 
are matrix elements for the SU(3) octet axial vector currents, mN = 0.939 GeV is the 
nucleon mass, yielding go = —0.21 GeV//a and g^ = — 0.17GeV//a. 



2.2. Expected axion flux from Li* -^"^ Li+a 

The decay of the flrst excited state of ^Li 

^Li* Li + 7 (0.478 MeV) (7) 

follows from ^Be electron capture (^Be+e~ — )■ ^Li* + z/g). Because this process can emit 
an axion of the same energy instead of a 7-ray and occurs for each neutrino emission, we 
can use the measured ^Be neutrino flux, to estimate the differential flux of axions 
arriving at Earth, as 



dE, 



Rq p 
J- - 



-exp 



2a{Tf 



(8) 



7 V27r(T(T) 

where Rq is the solar radius, (i$J^'^(r) is the ''Be neutrino flux at Earth emitted from a 
solar shell at radius r, k = 0.104 is the branching ratio of the ^Be electron capture to 
the flrst excited state of ^Li [33], cr(T) = E^^ykT/m is a thermal Doppler broadening 
of the emission line (about 0.2 keV at the solar core), and Fa/F^ is the axion-photon 
branching ratio. 

We use the values of 77 and (3 calculated in [30] as 77 = 0.5 and (3 = 1.0, thereby 
obtaining 

^ = 1.035((?o + ^73)'. (9) 

Since the resolution of the calorimeter in the energy region around 450 keV is a^ct ~ 100 
keV, we integrate over the Doppler broadening term and use the total neutrino flux at 
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Earth = 4.86 x 10^ cm ^ s ^ [M]. Thus, we wash out the Doppler term and obtain 
the total flux of 0.478 MeV solar axions 

<l>a = 5.23x 10^((7o + ^?3)' cm-V^ (10) 

It is interesting to note the relative insensitivity of $a to the choice of r] and 
P above. By varying 0.1 < rj < 0.9 and 0.6 < /3 < 1.4, $a varies between 
3.2 X 10^ < $a/(fi'o/3 + gsY < 8.8 X 10^ (in this range of r] and (3). Thus, the more 
critical parameter is the emission rate, which follows the ^Be neutrinos, $f 



2.3. Expected axion flux from D{p, '-f)^He 

Also of interest for hadronic axions is the radiative capture of protons on deuterium, 
referred to as proton-deuteron fusion [27]. The reaction 

p + d^3He + 7 (5.5 MeV) (11) 

occurs at a rate 0^ = 1.7x10"^* s~^ with only 1/3 of those being Ml transitions. In 
order to obtain the axion flux expected from this nuclear reaction we must evaluate the 
axion-photon branching ratio ([5]) using the correct values of parameters r] and /3, which 
is made difficult by the fact that f fTTj) is a 3-body nuclear decay. However, as pointed 
out by [27], ( fTTl) is predominantly isovector, implying that go/3 is very small and can be 
neglected. Under this assumption, equation ([5]) becomes [35] 

T\~2^Va^J i + s^\J^) ■ ^^^^ 

Using 6=0, ka ^ = 5.5 MeV, and /is = 4.71 we have 

^ = 0.98^3^ (13) 



Combining equation ( \T3h and the rate of proton-deuteron fusion reactions characterized 
as the Ml nuclear transition, we obtain the total flux of 5.5 MeV solar axions 

where ^0 is the Earth to Sun distance. 

Using the data obtained from the 6 month run in 2004, we can thus evaluate the 
sensitivity of the CAST 7-ray calorimeter to these axion signals. 



3. Data 



3.1. CAST and the high energy '-/-ray calorimeter 

The CERN Axion Solar Telescope (CAST) utilizes a helioscope design which exploits the 
increased axion-to-photon conversion probability for increased magnetic field strength 
and length as given by equation (j4|). The refurbished LHC dipole prototype magnet [36] 
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produces a nominal magnetic field of B = 9.0 T over a length of L = 9.26 m in each 
of the dipole's two 14.5 cm^ area magnet bores. The full system is mounted on a 
rotating platform with a vertical range of ±8° and an azimuthal range of ±40°. This 
range of motion allows for 1.5 hours of solar alignment during both sunrise and sunset 
year-round. The tracking system monitors the alignment of the magnet with the Sun, 
resulting in a pointing accuracy better than 0.01°. All remaining time is devoted to 
background measurements for the low-background X-ray detectors which are installed 
on both ends of the magnet. Until 2007, a conventional Time Projection Chamber 
(TPC) covered both magnet bores at one end to detect photons originating from axions 
during the tracking of the Sun at sunset. It was then replaced by two MICROMEGAS 
detectors, each attached to one bore. On the other side of the magnet, there is another 
MICROMEGAS detector covering one bore, and an X-ray mirror telescope with a pn- 
CCD chip as the focal plane detector at the other bore, both intended to detect photons 
produced from axions during the sunrise solar tracking. More details about the CAST 
experiment and detectors can be found in[I7l[T8l[32l[36l[371[3Hl[39lll0]. 

To cover a wide range of potential axion masses, the operation of the CAST 
experiment is divided into two phases. During the Phase I (2003-2004) [TTl [37] the 
experiment operated with vacuum inside the magnet bores and the sensitivity was 
essentially limited to ma < 0.02 eV due to the coherence condition. In the second 
phase (so-called Phase II) which started in 2005, the magnet bores are filled with a 
buffer gas in order to extend the sensitivity to higher axion masses. In the first part of 
this phase (2005-2006) '^He was used as a buffer gas. By increasing the gas pressure in 
appropriate steps, axion masses up to ~0.4 eV were scanned and the results of these 
measurements supersede all previous experimental limits on the axion-photon coupling 
constant in this mass range [18]. To explore axion masses above 0.4 eV, '^He has to be 
used because it has a higher vapor pressure than ^He. This allows for a further increase 
in gas pressure in the magnet bores to reach axion masses up to about 1 eV in the 
ongoing second part of Phase II that started in 2007 and is planned to finish by the end 
of 2010. 

Axions emitted in a particular nuclear decay channel will be essentially mono- 
energetic compared to the Primakoff spectrum expected from plasma processes. 
However, several candidate processes exist and so the corresponding mono-energetic 
axion lines are expected to be in the range from tens of keV to many MeV. The expected 
axion signal from any single channel is thus a coUimated beam of similarly mono- 
energetic 7-rays (a "peak") from the magnet bore during periods of solar alignment. 

The calorimeter consists of a cadmium tungstate scintillating crystal (CdW04 or 
CWO), which is also the type used in neutrinoless double-beta decay searches |41j . 
CWO offers good stopping power for 7-ray photons, very low internal radioactivity, 
good energy resolution and excellent pulse shape discrimination characteristics (see 
section 13.21) . The crystal is optically coupled to a light guide and photomultiplier 
tube (PMT) which is placed inside a lead-shielded cylindrical brass tube. This 
"tunnel" design maximizes signal acceptance and background rejection while respecting 
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Figure 2. (a) Schematic of tlic CAST magnet platform depicting the calorimeter 



placement behind the MICROMEGAS detector and the space constraints limiting 
the allowed shielding. |(b)| Schematic of the CAST high energy 7-ray calorimeter. 
The "tunnel" design with CWO crystal and PMT surrounded by modest amount of 
ancient and common lead passive shield. An active muon veto with dual PMT read-out, 
boratcd neutron absorber and continuous N2 radon purging complement the passive 
shielding. 



the space and weight hmitations on the CAST detector platform (see figure E^b)). 
These constraints also limit the allowed thickness (2.5 cm) of ancient and common 
lead shielding, which results in an elevated environmental background component 
compared to that achievable with fewer constraints. An active scintillating plastic muon 
veto, environmental radon purging with constant N2 fiow, a borated thermal neutron 
absorber, and a low-background PMT complement the minimalist passive shielding 
design. 

The large dynamic range and high stopping power for photons are necessary to 
achieve a good efficiency at high energies for a generic search. Detector components, 
shielding materials, and data processing were all designed in order to reduce the 
environmental backgrounds. Pulse shape discrimination (PSD) further reduces noise 
and events due to internal radioactive contaminations in the crystal. Finally, an LED 
pulser provides livetime monitoring. These square pulses are recorded and subsequently 
removed prior to analysis. 

Although similar searches have been conducted in beam dump experiments ^21 H3] , 
accelerators and terrestrial nuclear processes [H], CAST is the first high-energy- 
axion search using a helioscope. The order of magnitude increase in axion-to-photon 
conversion probability over previous helioscope searches and the increases intensity of 
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axion emission from the Sun as compared to accelerator searches makes the CAST 
calorimeter a very sensitive probe of low-mass pseudoscalars. Because axions serve 
as merely one example of such particles, a high-energy search should not be limited 
to only axions but should consider anomalous excess during solar tracking events 
generally [T21I21]. 

The data and results presented in this paper were obtained from CAST Phase I. 



3.2. Pulse shape discrimination and particle identification 

Following a method similar to [H], we have developed and applied a pulse shape 
discrimination (PSD) algorithm which exploits the distinct pulse shape characteristics 
of the CWO crystal in response to incident particle type. This algorithm relies on 
the difference in response for nuclear recoils and minimum ionizing particles due to 
the different mechanisms of scintillation for these two types of excitations. The long 
~20/is decay of CWO results from the presence of 2-4 separate decay constants [H] and 
enhances these differences and allows for efficient PSD. 

Particle calibrations were performed using alpha's (a's), neutrons (n's), and 
photons (7's) to create a weighting factor using a statistical average over many pulse 
shapes for each calibration source. The radioactive sources used for the calibrations 
represent a spectrum of energies: ^^^Am (a's), Am/Be and 252Qf ^^^^ ^^^^^ 7's), and 
multiple 7 sources (0.511, 0.662, 0.835, 1.173, 1.333, 1.836 MeV). Particle "templates" 
are formed with these data and are used to create a spectral weighting factor. The PSD 
algorithm is then applied to the data to remove backgrounds and noise. 

We construct the 7 template with 9322 events using the (1.8 MeV) source, the 
a template with 7781 events from the ^'^^Am source, and the n template with 450 events 
from the Am/Be source, each of which are shown in figure [31 Due to 7 contamination 
from the Am/Be source, several tests were performed in order to ensure a high purity 
of true nuclear recoil events. Layers of polyethylene shielding were added in steps 
to successively reduce the presence of neutron interactions in the crystal and isolate 
the 7 events in the sample. The pulse shape identification (PID) spectrum (shown 



in figure 3(b) for a and 7 calibrations and described in more detail below) was then 
remeasured in order to identify the region populated by the neutron recoils. Events 
were selected in this range for the measured neutron template. It is important to note 
here that these templates are obtained from raw data without signal filtering or artificial 
pulse shaping and represent the CWO response to mechanisms of energy deposited. 

The spectral weighting factor is determined using the particle templates described 
above and is defined as 



faih) - f,{tk) 



(15) 



fM + f\it,) 

Here, and are the pulse amplitudes as a function of time (tfc) for the a and 7 
templates, respectively. In order to select events based on PID for the n template, 
the PID distribution is first calculated using the weighting factor with the a particle 
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Figure 3. (a) Pulse shape "templates" for a's, n's, and 7's. The emphasis on the 



fast component of the peak is evident for the a and n recoils. The inset shows the 
peak region of each pulse shape (0-2 fis) and highlights the faster rise time observed 
for the neutron events as compared to a particles. |(b)| Pulsc shape identification (FID) 
separation for independent a and 7 particle calibrations. 



template as defined in equation [T51 Once the n template is formed using the selected 
events, may be replaced with f^. We refer to the factor in f|T5l) as a "spectral" 
weighting factor because it consists of a spectral shape which is to be convoluted with 
the entire waveform of each event during analysis. The full PID algorithm can be written 

as 



k=0 



fQ(tk) — f7(tk) 
fa(tk) + f7(tk) 



fevent(tk) (16) 



where we have converted this spectral shape in time into a single number, the event PID, 
by taking the sum over the kf time bins of the event weighted by the weighting factor. 
Every event is thus assigned a PID value based on this algorithm. To determine the 
optimum time window (kf) over which to integrate the pulse for a and 7 discrimination, 
we used the quantity 



(PID)„ - (PID)^ 



(17) 



which measures the peak separation of the PID distributions, where and are 
the RMS of the PID distributions measured during calibration. This quantity is then 
maximized for the PID distributions to obtain an optimum time window over which to 
integrate the pulse of tonsct < ^fe < 14/is, where tonsct is the rising edge of the CWO pulse. 
This window can be understood qualitatively given the calibration pulse templates in 
figure 3(a) which converge at approximately 10 fis. 

By using both the above described particle {a and 7) calibrations, event selection 
criteria were determined prior to the analysis of the signal data samples and remain 
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consistent throughout the analysis. These criteria are first set using 7 cahbration data to 
maintain a cahbration signal acceptance of greater than 99.7%, while rejecting electronic 
noise and square pulses from the livetime pulser. This yields a background a rejection of 
~50%, determined directly from the calibration data (figure 3(b) ). Due to the difficulties 
in obtaining the optimal algorithm for neutron rejection and the iterative procedure used 
for extracting the neutron pulse shape template, the neutron rejection capability has 
not been fully quantified for this analysis. The exact fraction of neutron recoils rejected 
would be exactly characterized with the help of a pure, mono-energetic neutron emitter 
which was not available at the time of the calorimeter commissioning and operation. 



3.3. Solar tracking and background data 

Both background and tracking events are considered for analysis using the same data 
quality criteria, while solar tracking events have the further requirement that the 
magnet be sufficiently aligned with the solar core. Corrections are then applied for 
a small background energy spectrum dependence on the pointing position of the CAST 
magnet, which is due to differences in natural radioactive background throughout the 
CAST experimental hall. By dividing the horizontal-vertical plane traversed by the 
CAST magnet into a set of cells, each of which represents a section of the wall or floor 
towards which the magnet points at a given time, the position dependence of all detector 
parameters is directly measured. 

In order to ensure the compatibility of the the final tracking and background 
spectra, we select background data only from those positions in which tracking data 
has been recorded and weight those data according to 

— m — 

where {dN^'^^'^ / dE)i is the background energy spectrum in the i^^ cell, is the 
tracking exposure time in the i^^ cell and {dN^'^^'^ / dE)^^ is the effective background 
after position normalization. Following these corrections, the background and tracking 
(signal) data sets can be reliably compared. 

The dataset includes a total of 1257 hours of total exposure time with 60.2 hours of 
solar alignment and 898 hours of background data. A summary of the statistics for this 
data set is shown in table [H The effective background data set following the position 
normalization procedure described above still consists of more than twice the tracking 
data, thus maintaining good statistics for background subtraction. 

To facilitate the analysis over such the 7-ray calorimeter large dynamic range in 
photon energies, the data are divided into three energy regions (0.2-3.0 MeV, 3.0- 
10 MeV, 10-100 MeV) and binned according to the detector resolution in each region. 
The energy spectra for both tracking and background in each energy range are shown 
in figures HI Environmental 7 radioactivity is very evident in the low energy region. 
7's from ambient activity (1.460 MeV) and ^osjl (2.614 MeV) from the ^ssjh 
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Table 1. Statistics of the data from CAST Phase I used in this analysis. 

Total exposure time 1257.06 h 

Solar tracking 60.256 h 

Background 897.835 h 

BCKG rate after cuts 1.429 Hz 

BCKG flux (above 200 keV) 0.1 cm^^ g-i 



e^e" annihilation (0.51 1 IVleV) 

'°K (1.460 IVIeV) 



'Bi (1.765 MeV) 




1 1.5 2 

Energy [MeV] 



Solar Tracking 

Background 




4 5 6 7 

Energy [MeV] 



(a) 0.2-3.0 MeV 



(b) 3.0-10 MeV 




20 30 40 50 

Energy [MeV] 



(c) 10-70 MeV 



Figure 4. Tracking and background energy spectra for each energy region in the 
analysis. The background spectra (narrow band) used have been normalized to the 
time and position of the tracking data (wider band) and PSD and energy cuts have been 
applied. Prominent in the low energy region are decays from ^"^K activity (1.460 MeV), 
208T1 (2.614 MeV) from the ^^^Th decay chain, ^"Bi (1.76 MeV) from the ^s^U decay 
chain, and e+e^ annihilation 7's at 0.511 MeV. 



decay chain exhibit prominent peaks in the data, along with e~^e annihilation 7's at 
0.511 MeV. 



4. Data analysis and results 
4.1. Expected axion signal 

Direct background subtraction from the tracking data permits the search for excess 
events in the residual energy spectrum. The expected signal from axion-photon 
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conversion is a coUimated "beam" of mono-energetic photons from the magnet bore 
during solar ahgnment. This results in Gaussian energy depositions in the CWO crystal 
for low energy (below 1.022 MeV) photons. 

Above 1.022 MeV, an axion conversion photon may pair-produce within the crystal. 
For each pair-production, there is the possibility that one or both annihilation photons 
escape. These annihilation escape peaks will lie at 0.511 and 1.022 MeV below the full 
energy peak and the efficiency for catching these events is characteristic of both the 
crystal and the energy of the incident axion-conversion photon. 

A standard MCNP4b [15] simulation of this spectrum for a 5.5 MeV photon, 
convolved with the detector resolution, is used to determine the calorimeter sensitivity 
to photons at this energy. Photon detection efficiency is nearly 48% when considering 
the entire range at 5.5 MeV for this signal (see table [2]). To validate these data, a 
laboratory replica of the MICROMEGAS X-ray detector which sits directly in front 
of the calorimeter in the experiment was constructed and the transmission efficiency 
through the detector material for photons of various energies was measured and 
compared to the Monte Carlo predictions. The data and the simulation were found 
to be in good agreement and the simulation was then used only to determine the peak 
efficiency and the relative peak heights for the deposition signal. 

The multi-peak signal shape and increased photon detection efficiency improves 
the sensitivity to excess events above 4.0 MeV. A general search along the entire energy 
spectrum of the calorimeter would require a full Monte Carlo analysis of the signal shape 
and its energy dependence. Here, only a 5.5 MeV photon signal has been investigated 
using this approach, while at all other energies below 10 MeV only a single Gaussian 
signal (corresponding to a full energy peak) is used. The 5.5 MeV signal is fit as two 
Gaussians, to a good approximation, which have a fixed peak-height ratio given by the 
simulation. The search for this signal is described in section 14.41 and the resulting fit to 
the data is shown in figure [71 

Above 10 MeV, photonuclear dissociation is both energetically possible and very 
probable, with cross sections near 1 barn for the tungsten and cadmium in the CWO 
crystal. This large cross section for interaction results in a much different signal shape 
than for low energy photons and can no longer be approximated by a Gaussian, which 
is taken into account for the six energies evaluated in the following section. In this 
energy regime, the total energy deposition efficiency and the signal shape is determined 
from a standard MCNP4b simulation and depends on the cross-section for photonuclear 
interactions above 10 MeV. Only 6 points above 10 MeV are evaluated in the model 
independent scan. 

4-2. Search for anomalous mono-energetic peaks 

We perform a generic search for excess photons across the entire dynamic range of 
the detector by fitting the known a mono-energetic signal shape to the residual energy 
spectrum remaining after background subtraction. For the analysis presented here, a 
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Figure 5. 95% CL limits on tlie number of counts above background yield limits on 
the flux of axion-conversion photons incident on the CAST calorimeter after correcting 
for detection efficiency and livetime. These limits may then be directly translated into 
the allowed flux of axions incident on the helioscope. Below 10 MeV, a Gaussian 
energy deposition signal is used. At higher energies, the photonuclear cross section 
alters the energy deposition signal and the analysis is performed for 20, 30, 40, 50, 60 
and 70 MeV. 

Gaussian energy deposition has been used below 10 MeV with the exception of the 
5.5 MeV photon signal, as stated in section HTTl 

The results of the search and extraction of 95% CL upper limits on excess photon 
flux are shown in figure [5l The structure present in the plot is a general consequence 
of statistical fluctuations in the residual spectrum of the data which lead to large or 
small 95% CL bounds on the Gaussian and are physically meaningful as they can point 
to incomplete subtraction or to slight statistical excesses in the data. Above 10 MeV, 
6 points were chosen at which to evaluate the presence of the photonuclear interaction 
signal (10, 20, 30, 40, 50 and 60 MeV). Although these data are difficult to interpret 
outside of the context of a particular model, they serve as a benchmark sensitivity for a 
relatively model independent hadronic axion searches using a helioscope, assuming only 
that the axion emission is mono-energetic. 

4-3. Calorimeter sensitivity to Li*-^'^ Li+a 

In order to evaluate the detector sensitivity to axion emission from specific decay 
channels, several parameters are included and are found in table O From MCNP 
Monte Carlo simulations, at the energy of the ^Li decay, 478 keV, the efficiency for 
peak energy deposition after accounting for the photon transmission efficiency through 
the MICROMEGAS detector is 56.8%. Including all known detector inefficiency, i^totah 
we can estimate the sensitivity to axions from '''Li decays using 

'^'478kcV = -Pa->7(fl'a7, nT-a., Eg., B, L) ^aTdct^total 
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Table 2. Summary of the calorimeter characteristics and data selection efficiencies 
used to calculate the expected number of photons incident on the detector for both 
''Li and D(p, 7)^He. 



Production channel ^Li D(p,7)'^He 



Peak efficiency 


^pcak 


56.8% 


47.5% 


Photon transm. efficiency 


^transm 


Energy resolution 


fdot 


99 keV (21%) 327 keV (6%) 


Livetimc 


Tdet 


93% 




Software cuts efficiency 


^cuts 


99% 




Conversion probability 


^ a— >7 


(see figure 


El) 




m, [eV] 



Figure 6. Axion-photon conversion coefficient P!^^^ versus axion mass. Since axions 
from nuclear processes have higher energies than Primakoff axions, the sensitivity for 
axion masses in case of ^Li decay (478 keV) and p-d fusion (5.5 MeV) is increased for 
1 and 2 orders of magnitude, respectively. 



= 1.368 X 10«P^^^(ma, E,, B, L)g%{go + g^f (19) 

Here we have factorized the conversion probabihty Pa^7 in equation (j4]) into the axion-7 
couphng constant g^^ and the remaining numerical term, Pj[^^(ma, -Ea, -B, L), which only 
depends on the magnet parameters and the axion energy and mass. Pa^^(ma, E.^, B, L) 
versus ma is shown in figure [6l 

To extract the signal, the background spectrum was subtracted from the spectrum 
obtained during solar alignment, and the resulting residual is shown in figure 7(a) Since 
no evidence of an axion signal was observed only an upper limit could be determined. A 
fit to the expected Gaussian signal shape yields a 95% CL upper limit on excess photon 
events at 478 keV of $478kev(95%CL) < 3.4 x 10"^ cm'^ s-\ By solving for g^^ in ^ 
we can write 



1_ / 'l'478kcv(95%CL) 1 

^^'-90 + 93^ 1.368 x 108 ^''^ 
The limits as a function of ma are shown in figure [8]for two different values of the nuclear 




Figure 7. Fits to tlie known signal shape after baekground subtraetion. (a) A single 
Gaussian fit for ^Li*^^Li+a results in a 95% CL upper limit for a 478 keV axion- 
conversion photon of $478kov(95%CL) < 3.4 x 10"'' cm-^ s-\ [(b)]A double Gaussian 
fit to the expected signal from D(p, 7)'^He yields the 95% CL upper limit for a 5.5 MeV 
axion-conversion photon of <I'5.5Mcv(95%CL) < 2.14 x 10^'"' cm~'^ s~^. 



coupling constant (fi'o + fi'a)- These values are chosen and evaluated using equation as 
representative of the range of couplings for a Peccei-Quinn scale of /a = 10^ — 10^ GeV. 

4.4- Calorimeter sensitivity to D{p,'~f)^He 

Using again the information from table [21 we can evaluate the limiting expression for 
(7a7 in the D{p,'j)^B.e channel using 

$5.5McV = 0.8878 X 10i°P^^^(ma, E,, B, L)g%gl cm-h-\ (21) 



The resulting spectrum after background subtraction is shown in figure 7(b) , For this 
reaction, the expected signal differs from that of ''Li since a 5.5 MeV 7-ray can pair- 
produce within the calorimeter. This escape peak structure has been taken into account, 
including the fixed peak-height ratio, resulting in a 95% CL limit on excess photons of 
$5.5Mev(95%CL) < 2.14 X 10"^ events cm"^ s~^. Using equation (12T1) we have 




5.5Mev(95%CL) 1 

0.8878 x 1010 P^:^^' ^ ^ 

We again use the two different values of the nuclear coupling constant (g^), 
evaluated using equation and a Peccei-Quinn axion scale of /a = 10^ — 10^ GeV. By 
plotting the axion-photon coupling constant versus axion mass, we see that the limits 
are weaker than those obtained by the CAST X-ray detectors in 2003. 

5. Conclusions 



The CAST photon calorimeter provides a search for high-energy axion-photon 
conversions during periods of solar alignment. This is the first such search for high- 
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ma [eV] 

Figure 8. The limits obtained on the axion-photon couphng versus axion rest mass 
for 478 keV axions from ^Li decay and 5.5 MeV axions from proton-deuteron fusion 
for two values of the nuclear couplings. This parasitic axion search has not found any 
evidence for new pseudoscalar particles coupling to nucleons. 

energy pseudoscalar bosons with couplings to nucleons conducted using a helioscope 
approach and provides an important cross-check for other searches focused on nuclear 
decay, such as [301 US]- Furthermore, as discussed in [27], the search for pseudo-scalar 
emission from proton-deuteron fusion (D(p, 7)^He) is potentially sensitive to a more 
general class of new particles than only Primakoff or hadronic axions due to the presence 
of both Ml and El transitions and can couple particles of various spin-parity. 

In making use of the CAST magnet for an axion search strategy not initially 
foreseen, the achievable sensitivity is severely limited and a dedicated high-energy axion 
search performed underground and without shielding limitations would be able to reach 
background levels many orders of magnitude lower. Such levels are necessary to reach 
the very small axion flux expected from the two axion emission channels considered in 
this search. CAST remains a unique instrument with unprecedented sensitivity allowing 
for new searches for anomalous solar emissions in the form of new axion-like particles 
with coupling to photons. 
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